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Abstract 

We study observational implications of the stochastic gravitational wave back- 
ground and a non-Gaussian feature of scalar perturbations on the curvaton mech- 
anism of the generation of density /curvature fluctuations, and show that they can 
determine the properties of the curvaton in a complementary manner to each other. 
Therefore even if Planck could not detect any non-Gaussianity, future space-based 
laser interferometers such as DECIGO or BBO could practically exhaust its param- 
eter space. 



1 Introduction 



The idea of inflation [1] has become a standard paradigm since it was proposed in the 
early 1980's. Inflation is an accelerated expansion epoch in the very early Universe, which 
solves the flatness, the horizon and also the monopole problems naturally. Furthermore 
quantum fluctuation of the inflaton 0, which is a scalar field responsible for the accelerated 
expansion, can provide the seed of cosmic density/curvature fluctuations [2] observed 
through the cosmic microwave background (CMB) anisotropy [3], galaxy clustering, etc. 

Another prediction of inflation, which is in fact more generic, is the generation of the 
stochastic gravitational wave background or the tensor perturbation, whose spectrum is 
nearly scale-invariant [4-34] (see Ref. [35] for a review). The amplitude of the gravitational 
wave is simply proportional to the Hubble parameter, H ini , during inflation and many 
inflation models predict detectable amplitude of gravitational waves. In particular, future 
space-based gravitational wave detectors such as DECIGO [36] and/or BBO have a chance 
to detect inflationary gravitational wave background. Recently, it was pointed out that 
thermal history of the early universe is imprinted in the spectral shape of the gravitational 
wave background [14,21,28,29,31,32]. In particular, the reheating temperature of the 
Universe after inflation can be determined or constrained from future gravitational wave 
experiments [29,31,32]. Thus any detection of primordial gravitational wave background 
gives useful information on the early Universe. 

Traditionally, the spectrum of the comoving curvature perturbation is parametrized 

as 



Here k* is the pivot scale which we take k* = 0.002 Mpc -1 and n s is the scalar spectral 
index. Inflation predicts nearly scale-invariant spectrum with n s = 1, with its precise 
value determined by the shape of the potential [2], and it agrees well with observations so 
far [3]. If (j) has a single component with a canonical kinetic term, the resultant curvature 
fluctuation is Gaussian distributed, which is again in agreement with the observation 
today [3]. 

However, the origin of the density perturbation is not limited to the quantum fluc- 
tuation of the inflaton. Another scalar field, called curvaton, may be responsible for the 
generation of the observed density perturbation [37,38]. It is undoubtedly an important 
task to distinguish them in order to understand the physics of the early Universe. One 
such signature may come from the non-Gaussianity in the CMB anisotropy [39] , since the 
curvaton scenario can produce large enough non-Gaussian feature to be detected [40,41], 
while standard inflation models predict negligible non-Gaussianity [42-45]. However, the 
situation where the curvaton generates large non-Gaussianity is somewhat limited and 
it is possible that the curvaton accounts for the observed density perturbation without 
generating large non-Gaussianity. 

In this paper, we point out that the observation of stochastic gravitational wave back- 
ground plays a very important role to probe the physics of the curvaton scenario. In the 
case that the curvaton generates negligible non-Gaussianity, there must be an entropy 
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production process by the curvaton decay itself, and such a non-standard thermal history 
is imprinted in the spectrum of the inflationary gravitational wave background. Future 
space-based gravitational wave detectors, such as DECIGO and/or BBO may be able to 
do this job. This opens up a possibility to find an evidence of the curvaton scenario. 
Thus measurement of the non-linearity parameter /nl which is a simplified measure of 
non-Gaussianity and that of tensor perturbations by DECIGO/BBO play complementary 
roles to each other. 

The rest of the paper is organized as follows. In Sec. [21 the gravitational wave back- 
ground spectrum and non-Gaussianity in the curvaton scenario are summarized. In Sec. [3] 
the detection possibility is discussed. Sec. 0] is devoted to the conclusion. 

2 Features of the Curvaton Scenario 

In this section, we summarize the cosmological consequences of the curvaton scenario, 
namely, the spectra of gravitational wave background and the density/curvature pertur- 
bation, and non-Gaussianity. All these ingredients are essential for probing the curvaton 
scenario and constraining the parameter space, as will be discussed in Sec. 131 

2.1 Gravitational wave background spectrum in the curvaton 
scenario 

In the inflationary era, quantum fluctuations of the RMS amplitude of H in f/(27r) are 
induced on all massless fields in each Hubble time. Here "massless" means that the mass 
is much smaller than H-^. The tensor perturbation of the metric consists of two free 
massless scalar components, which can be quantized in de-Sitter space-time. After a 
mode left the horizon during inflation, it actually becomes a classical fluctuation which 
can be viewed as a stochastic gravitational wave. Thus inflation necessarily generates 
a gravitational wave background from a cosmological scale to a sub-kilometer scale The 
former is a target of detection by the B-mode polarization of CMB anisotropy and direct 
detection by space-based laser interferometer experiments. 

Tensor perturbation of the metric is defined by the following line element, 

ds 2 = a(r) 2 [-dr 2 + (<% + h^dx^x^ , (2) 

where a(r) is the scale factor and hij denotes the metric perturbation satisfying the 
transverse-traceless conditions d l hij = and h\ = 0. Thus has two physical degrees 
of freedom, which are denoted as h x with A = +, x. In the inflationary era, the tensor 
perturbation has a quantum fluctuation whose spectrum is given by 
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where the tensor spectral index is given by n t = — 2e. Here e is one of the slow-roll 
parameters during inflation defined by 



VJ ' ^ = M ^ (4) 

where V is the inflaton potential and the prime denotes derivative with respective to the 
inflaton field 0, and Mp\ = (87rG) -1 / 2 is the reduced Planck scale. 

After the production of the gravitational waves during inflation, the amplitude of each 
Fourier mode remains constant when the corresponding mode lies outside the Hubble 
radius. However, once it enters the horizon, its amplitude decreases as oc a -1 . Thus the 
present energy density of the gravitational wave background per logarithmic frequency 
interval is written as 

d Pgw k 2 A ( P )2 m ( a-m{k) 



dink 327rGag Afc { a ) ' (5) 

where ao is the present scale factor, and a- m (k) denotes the scale factor at which the 
corresponding mode with wave number k enters the horizon. It behaves as oc k~ 2 (k°) 
for the mode which enters the horizon at matter (radiation) dominated era. Therefore 
thermal history of the Universe is imprinted in the spectrum of the gravitational wave 
background, and this is the reason why the observations of the gravitational wave are 
expected to have great impacts on cosmology [14,21,28,31,32]. 
In terms of the density parameter, it can be rewritten as 

^(k) = 7^2^)' ( 6 ) 



where H is the present Hubble parameter, and 



4/3 



A£(A0 = A« a (*R ( 3 JlM) f 9 -^) ( ) 2?W2?W, (7) 

V z k J \ 9*0 J \g*s{Ti a )J 

where g*(T in ) denotes the effective relativistic degrees of freedom at the temperature T in 
when k-mode enters the horizon, and ji(z) is the spherical Bessel function of the first 
rank with Zk = 2k/(aoH )^ The transfer functions Ti(x) and T 2 (x) are given by [8,32] 

T 2 (x) = 1 + 1. 57x + 3.42x 2 , (8) 
T|(x) = [1 - 0.32x + 0.99X 2 ]- 1 . (9) 

The former connects the gravitational wave spectrum of the mode entering the horizon 
before (x eq = k/k eq > 1) and after (x e „ < 1) the matter-radiation equality, where k eq = 
a(t eq )H(t cq ) = 7.3 x 10 _2 fi m /i 2 Mpc _1 o The latter transfer function connects the mode 



1 There was an error in Eq. (14) of Ref. [32]. Eq. ([7]) is the correct one. 

2 The effect of neutrino free streaming is known to lead a suppression on the gravitational wave 
background spectrum around the frequency ~ 10~ 9 Hz [15,26]. But our concern is around 1 Hz, and 
hence we simply neglect this effect. 
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entering the horizon before (x R = k/k R > 1) and after (xr < 1) the reheating subsequent 
to inflation. Here k R is the comoving wavenumber corresponding to the horizon scale 
at the reheating epoch when the Universe became radiation dominant. Without any 
significant entropy production after the reheating, it is given by 



A:r -1.7x10 13 M P c- 1 ^* s(Tr) V / V T * 



o V 106 - 75 / V 106 GeV 

10 ^ 1/4 



(10) 



with Tr and being the reheating temperature after inflation and the decay rate of the 
inflaton, respectively. This corresponds to the frequency / =0.026Hz for Tr = 10 6 GeV, 
which is close to the most sensitive frequency range of the planned future space-based laser 
interferometer experiments, DECIGO or BBO. Thus by observing the spectral shape of 
the gravitational wave background, the reheating temperature of the Universe can be 
determined or constrained. 

As mentioned above, the above correspondence assumes the standard thermal history 
with no significant entropy production after reheating. In the curvaton scenario, however, 
its coherent oscillation may once dominate the Universe, which introduces an additional 
matter-dominated era, and then decays releasing huge amount of entropy. In this case 
the gravitational wave spectrum receives an additional suppression [14,31,32] which can 
be quantified by the dilution factor, F, defined by 



a? Xi 



a (T>3 (rg) _ J ^kf for m„>T, 
s(T R )a 3 (T R ) 1 T osc > {lL) 

where is the decay rate of the inflaton, T a is the radiation temperature just after the 
curvaton decay, and T osc is the temperature at which the curvaton begins to oscillate. 
This expression is valid for F>1. The resultant gravitational wave spectrum is given by 

Al(k)=A^\k)nl( 3jl{Zk) ^ ( 9 * {Ti ^ ( 9 * s0 Nl '° 



Zk J \ g*o J \g*s(T in ) J (12) 

x Tl (x eq )T 2 2 (x a )T? (x aR )T% (x R (F) ) , 

where x a = k/k a with k a given in an analogous way with (ITU]) after replacing Tr with T a , 
and x aR = k/k aR with k aR = k a F 2 / 3 and x R (F) = k/k R (F) with k R (F) = ^F' 1 / 3 . 

In Fig. [1] we show spectra of the gravitational wave background as a function of its 
present frequency. In the top panel, the spectra for if in f = 10 14 GeV and 10 13 GeV with 
the reheating temperature Tr = 10 6 GeV are shown. Also plotted are sensitivities of 
DECIGO with a correlation analysis (blue dashed line), ultimate-DECIGO (red dotted 
line), and correlation of analysis of ultimate-DECIGO (purple dot-dashed line) [20]. In 
the bottom panel, the gravitational wave background spectra in the presence of entropy 
production is shown, for F — 10 and T CT =10 GeV and Tr = 10 7 GeV. 
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Figure 1: (Top) Spectra of the gravitational wave background for inflationary scale H m f = 
10 14 GeV and 10 13 GeV. Here we have taken Tr = 10 7 GeV. Also shown are sensitivities 
of planned space-based gravitational wave detectors, DECIGO with a correlation analysis 
(blue dashed line), ultimate-DECIGO (purple dotted line), and correlation of analysis of 
ultimate-DECIGO (red dot-dashed line). (Bottom) Same as the top panel for the dilution 
factor F = 10 for 7^=10 GeV and T R = 10 7 GeV. 
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2.2 Scalar perturbation in the curvaton scenario 

The curvaton is a scalar field other than the inflaton, which remains light during inflation 
and has a quantum fluctuation 5a ~ H- m i/(2ir). Then, according to the 8N- formalism [46], 
the comoving curvature perturbation 1Z C is given by 

K c = NJa+^N a(7 (8a) 2 , (13) 

up to the second order in 5a, where N is the local number of e-folds, given by the integral 
of the local expansion from an initial spatially flat hypersurface to a final uniform density 
hypersurface, and N a is its derivative with respective to a. Note that we assume the 
curvaton behaves as free scalar field and hence 5a can be regarded as a random Gaussian 
variable. Then we obtain [47] 

where cr, is the initial amplitude of the curvaton during inflation, and 

3p CT 



R 



4p r + 3p a 



(15) 

<r decay 



roughly denotes the fraction of the curvaton energy density to the total energy density at 
the epoch of curvaton decay. Here p a and p r are the energy densities of the curvaton and 
that of radiation, respectively. In order to reproduce the observed density perturbation 
of the Universe, we must have 

w=?S)= ,x ^ « i6) 

Notice that the inflaton also generates curvature perturbation whose magnitude is 
^n c 6 = ^taf / (87r 2 eM| 1 ) . It is much smaller than curvaton's contribution (1TB]) by assump- 
tion. This leads to a constraint 

9 ( Oi \ 2 _ . 3 / i^mf 



where we have used the WMAP normalization (fi~6l) . We can calculate the tensor-to-scalar 
ratio r in the curvaton scenario, as 

Al c (k.) IP\MpJ VlO u GeV^ ' [ 1 

It can be checked that this is smaller than 16e using Eq. ( TTTT) . which is the prediction of 
the standard inflation scenario, once we assume that the curvature perturbation from the 
inflaton should be smaller than that from the curvaton. 
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The scalar spectral index is given by [38] 

2m 2 

n s = l-2e + --f, (19) 

where m a is the curvaton mass. We can determine it through the relation 

ml = ^(n s -nt- l)Hlt = ll(n s - n t - l)r x (10 14 GeV) 2 , (20) 

in principle. But it is generically much smaller than H? nf , so that it would be more 
practical to use the above equality as a consistency relation, 

n.-nt = l, (21) 

for the curvaton scenario. 

Thus in the curvaton scenario, both the scalar and tensor spectral tilt are predicted 
to be — 2e, and hence a measurement of the tensor spectral index by future space-based 
gravitational wave detectors [19] will give an evidence of the curvaton scenario. 



2.3 Non-Gaussianity in the curvaton scenario 

Statistics of the observed CMB anisotropy is currently consistent with Gaussian distri- 
bution. The deviation from Gaussianity is parameterized by the non-linearity parameter 
/nl, whose definition is given by 

n c = + jj/ N L^ c (s)2 , (22) 

where 1Z C ^ denotes the Gaussian part of the curvature perturbation. WMAP5 result 
gives a constraint on it as —9 < /nl < HI at 95% C.L. [3]. From Eq. ( 1T4"1) . we can 
estimate the non-linearity parameter ag^| 

^ = 3S - > - H ' < 23 » 

Thus /nl can be significantly large for small R, and this may provide an observational 
hint of the curvaton scenario if large /nl is detected, because standard inflation models 
predict a small non-linearity parameter, /nl ~ 0(e,rj). Here we have assumed that the 
curvaton potential is quadratic, V = (l/2)m 2 r a 2 . If it deviates from the quadratic one, 
the prediction of /nl changes [41] and in this case the coefficient of the trispectrum, #nl, 
may be useful to distinguish curvaton models. We do not go into the detail on this point. 

Notice that in the case where a large non-Gaussianity (/nl ^ 10) is obtained, the 
curvaton decay does not increase entropy because it must be a subdominant component 

3 Here we assume that there are no CDM/baryonic isocurvature perturbations, which might cause 
other types of non-Gaussianity [48-50]. 
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at the instance of its decay (i.e., R <C 1), so that the gravitational wave background 
spectrum is not modified. In this sense, the detections of /nl by the CMB observation 
and entropy production process by the gravitational wave background as a probe of the 
curvaton scenario are complementary to each other. 

To summarize, we have five possible observable quantities: scalar spectral index n s , 
tensor-to-scalar ratio r, tensor spectral index n t , the dilution factor F and non-linearity 
parameter /nl- Among them, r will be accurately determined by future B-mode polariza- 
tion measurements. Direct detection of gravitational waves will observe n t , and comparing 
it with n s will confirm the curvaton scenario. Then, either large enough F or /nl will be 
observed depending on whether the curvaton once dominated the Universe or not. In the 
former case, F can be determined by comparing r and directly observed magnitude of the 
gravitational waves. In the latter case, /nl will be determined by CMB measurements 
such as Planck. Then these observables may be used to pin down a curvaton model. 

In the next section we investigate the possibility to detect either signature of the 
curvaton scenario and how they can fix properties of the curvaton. 

3 Probing the Curvaton Scenario 

We have seen that a curvaton scenario may leave distinct signatures on either a shape of 
the gravitational wave background or primordial non-Gaussianity. Here we show param- 
eter regions where a curvaton scenario has characteristic features on either of them. 

Curvaton models are characterized by two parameters, the initial amplitude of the 
curvaton o"j and its decay temperature T a . The relevant quantity is the abundance of 
the curvaton coherent oscillation at the time of its decay, which depends on these two 
parameters and the reheating temperature after inflation, Tr, but does not depend on 
the curvaton massfl Once these parameters are fixed, we can estimate the inflation 
scale H m t in order to reproduce the observed magnitude of the density perturbation (see 
Eq. (fl6l) ). The inflation scale H m { also gives overall normalization of the gravitational 
wave background spectrum, which can be directly measured independently using B-mode 
polarization of CMB. The spectral shape is determined by Tr and the curvaton abundance 
at its decay. Thus we can uniquely predict the gravitational wave signal at an arbitrary 
frequency for each parameter set (Tr, T CT , <jj). Similarly, the level of non-Gaussianity 
depends only on the curvaton abundance at its decay, as given in Eq. (1251) . and hence is 
also uniquely predicted. 

Figures 2 depict the parameter region of the curvaton which can be probed by obser- 
vation of DECIGO/BBO. The upper panel represents the region accesible with a single 
ultimate-DECIGO and B-mode experiments such as EPIC [51], CMBPol [52], and Lite- 
BIRD [53] which are expected to reach down to r ~ 10~ 3 . The lower panel shows the 
ideal case of correlation analysis of ultimate-DECIGO together with low-noise delensed 

4 If the curvaton begins to oscillate after the inflaton decays, the mass dependence appears. In this 
case, however, it is sufficient to replace 7r with T osc , the temperature at which the curvaton begins to 
oscillate, as in Eq. (fTTj) . 
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CMB map which would hopefully reach r~2x 10~ 6 [54] (see also [55]). 

Sensitivities of these direct detection experiments for the low frequency are limited by 
stochastic noise from white dwarf binaries [56], and hence we have cut the sensitivities 
below 0.1 Hz. We have set e to satisfy the constraint (fT7|) . but the precise value of e does 
not affect the results as long as it is sufficiently small. The region above the purple wedge 
corresponds to r > 0.2 which is excluded by WMAP. 

Also plotted there are contours of the non-linearity parameter /nl and the dilution 
factor F. The region with /nl > 100 is also disfavored by WMAP. On the other hand, 
Planck can measure it if it lies in the range 10 < /nl- As argued before, it occupies a 
detached domain from the region with F > 1. Therefore below we study how the curvaton 
parameters are determined in each domain separately. 

First let us consider the case /nl ^ 1 is confirmed by, say, Planck experiment. Then 
from (123]) we find R = 5/ (4/nl) and therefore (fl8|) determines the initial amplitude of 
the curvaton in terms of the observable quantities alone as 

" 5 < rVI \ (24) 



M P1 12/ NL \2 

In this case the curvaton decays before dominating the cosmic energy density. Hence we 
find 

i / \ 2 ,m \ i / \ 2 



P 



a decay 



6 \M VI ) a(T R ) 6 \MpiJ T r _ 



These two equalities lead to 

T a 5r 



If <jj \ a(T„) 1 / Oi \ Tr 4 . . 



(26) 



Tr 576/ nl 

Thus both of the curvaton parameters are fixed by the observable quantities in this case. 
Finally we note that in this case there is a simple relation between T ct /Tr and cii/Mpi as 



T a /nl / °. 



T R 10 VM P1 



(27) 



which was used to draw the contours of /nl in Figures 2. 

Next we consider the case the curvaton dominates the energy density of the universe 
when it decays releasing significant amount of entropy with R = 1 and F > 1. Then from 
PSD we find 

a,- / r \ 1/2 

(28) 



Mpi VI 
and from (TTTT) we obtain 

T r 

— = — ^— , (29) 
T R 108F' v ; 

for T osc > Tr. We can determine F if we can measure tensor perturbations both by the 
B-mode polarization of CMB and DECIGO/BBO, provided there is no other source of 
entropy production mechanism after reheating besides the curvaton. Then all the relevant 
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parameters are again fixed by the observable quantities. The corresponding parameter 
space is shown by the orange regions in Figs. [2J In case there is other entropy production 
besides curvaton, the right-hand-side of (1291) will give a lower bound on T ct /Tr. Finally, if 
B-mode polarization measurements fail to find the tensor mode and if it is detected only 
by DECIGO/BBO, which is shown by the yellow regions in Figs. [2j the values of F and r 
cannot be determined independently. Only the combination rF~ 4 ^ is determined in this 
case. 

4 Conclusion 

In this paper we have investigated a possible observable signatures from curvaton scenar- 
ios, including non-Gaussianity, detection of tensor modes in CMB and direct detection of 
gravitational wave background. First the measurement of the tensor power spectrum by 
CMB is very important to identify the curvaton scenario through the consistency relation 
(|2ip . If the curvaton once dominates the Universe, an entropy production process by its 
decay is imprinted in the gravitational wave background spectrum and can be confirmed 
by future space-based laser interferometer experiments. For the opposite case where the 
curvaton is a subdominant component before it decays, a large non-Gaussianity is pre- 
dicted, which may be confirmed by Planck experiment. In this sense, the gravitational 
wave signal and non-Gaussian signal are complementary to each other. This allows us to 
probe the curvaton scenario for large parameter spaces, giving information on the prop- 
erties of the curvaton such as its decay rate and amplitude, etc. It will in turn provide 
important information on high-energy physics beyond the reach of terrestrial experiments. 
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